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ABSTRACT 


In this work, the nucleon-nucleon interaction at 
low energies (below 300 MeV) is examined in a non-pertur- 
bation-theoretic manner. Two models are considered, the 
first being that where the interaction is mediated by scalar 
mesons (the scalar model), and the second, by pseudoscalar 
mesons (the pseudoscalar model). 

The scalar model is treated first. Equations for 
the scattering and bound states are derived in configuration 
Space. Explicitly non-local and energy-dependent potentials 
are obtained. In momentum space, an equation for the half- 
off-shell T-matrix is derived. The Born approximation for 
this model is then obtained. Half-off-shell T-matrix 
elements, and the phase shifts up to 230 MeV are calculated. 
For the bound state, a variational approach is taken to 
obtain a value of the coupling constant that gives the 
desired deuteron binding energy. 

The pseudoscalar model is accorded a similar treatment. 
There are four scattering channels in this case, and par- 
ticular attention is paid to the S = 0, T = 1 channel; it is 
treated analogously to the scattering state in scalar theory. 
The T-matrix is once again obtained. Finally, for the bound 
state, viz. the deuteron, configuration space potentials 
are obtained, and the equations for the s- and d-waves 


are derived. 
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I INTRODUCTION 


It had become apparent in the late 1920's, that 
nuclei could not contain electrons, because of their size, 
spin, statistics, etc. Following the discovery of the 
neutron in 1932, it was established that all nuclei were 
composed of neutrons and protons - known jointly as nucleons. 

The next step, obviously, was to try to understand 
the nature of the forces that keep a nucleus intact. As the 
simplest, and in many respects, the most fundamental of all 
nuclei, is that consisting of two nucleons, such systems 
have undergone a great deal of surveillance and study, both 
experimentally and theoretically. A great deal of data on 
the two-nucleon system has been amassed, at various energies, 
and large quantities of theoretical work have accumulated in 
the literature. One reason for this is that it is hoped 
that perhaps understanding N-N forces would lead to further 
clarification of, and more tractable methods of handling, 
multi-nucleon systems, and also nuclear structure problems. 

In this work, the 2-nucleon problem is taken up 
once again. The nucleons are non-relativistic, and hence 
the models will not be valid much higher than the single 
pion threshold (or about 300 MeV lab energy). In Section 
II, a very brief historical resumé is given of some of the 
many approaches taken. As this work involves a certain 
technique of quantisation, Section) Tilt is devoted to out- 


lining it; the significance of some of the terms that 
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appear will also be pointed out. Section IV is concerned 
with the application of this technique to the first of 
two models, viz. the model with scalar meson exchange. 
Both momentum and configuration space representations are 
discussed, and equations for the scattering amplitude and 
the T-matrix are derived. Certain physical quantities are 
then extracted (within an approximation). A similar treat- 
ment is accorded the theory with pseudoscalar mesons, in 
Section V, and phase shifts are extracted, again within the 
limits of the technigque of evaluation. Finally, the con- 
clusions and a discussion of results are contained in 
Section VI. 

The purpose of Appendices A - J is to remove the 
somewhat lengthy calculations (at least in some cases) 
from the main body of the text, thereby retaining some 


degree of continuity in it. 
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II HISTORICAL BACKGROUND 


The problem of the interaction between two nucleons 
is one of the most fundamental in all of physics. Over 
the past four decades, it has been dealt with by a large 
number of physicists in a multitude of ways. Whereas, in 
1960, the situation was such as to prompt the following 
quote Of MUL. Goldberger). 

"There are few problems in modern theoretical 

physics which have attracted more attention 

than that of trying to determine the funda- 

mental interaction between two nucleons. It 

is also true that scarcely ever has the world 

of physics owed so little to so many. In 

general, in surveying the field, one is 

oppressed by the unbelievable confusion and 

conflict that exists", 
the last fifteen years have to be considered as having 
been fairly productive. 

There are a number of reasons why the N-N inter- 
action warrants study; among them is the fact that for 
low energies, the N-N interaction is currently the most 
well-understood strong interaction process. Hence it has 
to be looked upon as a proving ground for ideas in ele- 
mentary particle physics ‘2? , Also, an understanding of 


N-N interaction would facilitate the treatment of nuclear 


structure calculations to a great degree. 
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The purpose of this section is to outline, in 
brief, the history of such work. It has, of necessity, 
to be brief, because of the very large number of articles 
in the literature relevant to the problem. A fuller 
review would be very voluminous, and apart from the fact 
that such reviews are published from time to time, is not 
of prime importance here. The purpose of this section is 
to try and view the N-N interaction in itsS proper pers- 
pective, and to describe, albeit briefly, the major con- 
tributions toward solution of the problem, their relative 
merits, and their shortcomings and pitfalls, if any. 

Basically, there have been two different approaches 
taken, viz. (a) phenomenological potential models, and 


(b) field-theoretical models. 


II.1 Phenomenological Models 


The prime motivation behind the use of potentials 
to describe what is, in fact, a strong interaction originally 
was the relative success of such an approach in electromag- 
netic and gravitational problems. It is phenomenological 
in that the potentials are constructed so as to reproduce 
existing data from experiments. 

One of the earliest attempts was by Gammel, Christian 
and Thaler ‘3) , and it was to fit all known data by a 


potential of the form 


V(r) = Vo (r) a Vp (5) S15 CLT etl) 
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The subscripts C and T stand, of course, for the central 
and tensor components of the potential, respectively, and 


Si9 is the usual tensor operator 
> “a > “A 
Si> = 3(0,-r) (05-4) - 0 


Eqn. (II.1) in fact represents the most general momentum 
independent potential which is local, conserves total 
angular momentum, parity and charge, is time-reversal in- 
variant, and charge independent. However, a good fit was 
not possible, and in order to improve the situation, the 
restriction of momentum independence was removed; a term 
linearly dependent on the relative momentum was added 


(4) 


(which leads to the Wigner type ‘of (potential) to give - 


V(r) = Va(r) + Viplr)S,5 + Vig (x) L-$ (Ta) 


Although with (I1I.2), fits were improved, they were still 
not satisfactory. In fact, if one assumes that the momentum 
dependence is purely linear, then only one potential V(r) 


can act in singlet-odd states, and another in all singlet- 


(4) 


even states . However, it can be shown that it is in 


fact impossible to describe all singlet-even states by a 


Single central potenti aier s 


Subsequent to attempts of this nature, have come 


three of the most widely used potentials. In 1962, the 


(6) (7) 


Hamada-Johnston and the Yale potentials were con- 
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structed; neither was of the Wigner type. Both of these 
were vastly improved, in terms of data fits, over previous 
attempts; the price that was paid, was, of course, sim- 


plicity of form. The potentials are: 


> > 

Vig (©) = Ve (+) + Vip (2) S45 + Vig (L)L-S + Vin (MP2 
with 
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The Hamada-Johnston potential, for example, reproduced 
low energy (<315 MeV) 2N data better than any other potential 
to date; data such as effective range expansion parameters, 
deuteron parameters etc. Vig (¥) is a hard-core potential, 
with its functional form outside the core being a combina- 
tion of Yukawa potentials. 

The third potential which has met with reasonable 


(5). 


success is the potential of Reid This too, “Ls 'non— 
local, in the sense of momentum dependence, and is of a 


different type in each angular momentum state. In addition, 
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both hard and soft cores can be made to fit the data; 
again, the radial dependence is of the Yukawa type. 
However, as with Ving and Vvale’ 


quires a fairly large number of parameters; this, of course, 


the Reid potential re- 


is an obvious consequence of attempts to fit different 
potentials in different angular momentum channels. 
An even more recent development is that of de-Tourreil 


(8) 


and Sprung Their potential is a "super-soft core" 


potential, and is of the following form: 


Like the Hamada-Johnston and Reid potentials, this one also 
approaches the OPEP (one-pion exchange potential) tail at 
large r. The three differ from each other near the origin - 
hence the names "hard-core", "soft-core", and "super-soft 
core". However, the fit to data up to 300 MeV is fairly 
good for this potential as well. 

Briefly, the current status of phenomenological po- 
tentials is the following: they all give (more or less) 
the same on-shell parameters such as phase shifts. And 
this is not sufficient to determine the validity of any 
particular potential. This is related to the fact, that 
apart from problems of parameter fitting etc., there is a 
more basic difficulty with regards to the phenomenological 


approach. Elastic processes are, in fact, insufficient to 
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determine, or test fully, an N-N potential ‘9, Even if 

one has a theoretical two-nucleon potential which fits 

all phase shifts (which, in practice, one will never be 

able to obtain experimentally), one has to consider, in 
addition, off-shell phenomena such as nuclear bremsstrahlung. 
This is because, given a potential which gives a set of 
phase shifts for all angular momenta and energies, it is 
possible to construct a high order of infinity of potentials 
which will reproduce this data i.e. the class of phase 
equivalent potentials. Phenomenologically constructed 
potentials will therefore have little to say about off- 
shell phenomena. However, the current experimental work 

has; nOtavyet been able ctoiconclusively test off-schell phe- 


nomena. 
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II.2 The Field-Theoretical Approach 


The basic philosophy behind this approach (in fact, 
behind any field theory) is that the interaction of the 
system may be described by the exchange of particles. Due 
to the relationship between the mass of exchanged quanta, 
and the range of the force mediated by this exchange, it is 
possible, in practice, to describe the interaction at 
varying distances from the centre of the system, by the 
exchange of particles with the appropriate mass. A useful 
guideline in this respect has been the approach of Taketani, 


Nakamura and See WEOne 


The region of interaction is divided 
into three ranges: the classical region (r 2 2f), the dy- 
HWamical region, (Zt oc. r >) jaf), ‘and? the icore i(r4<"0L)..: In 

the classical region, it is expected that one pion exchange 
(OPE), by virtue of the pion mass, will dominate. When one 
gets into the intermediate (or dynamical) region, exchanges 
of other types of particle begin to enter; finally, in the 
core region, the situation is complicated enough to warrant 
the use of a phenomenological core. 

The one-pion exchange potential (OPEP) tail is the 
only well-established meson-theoretic potential; in fact, 
its correctness has been established experimentally since as 
far back as the late fifties.,, Most, if not all, of recent 


work has gone into explaining the behaviour in the dynamical 


region, where one no longer has OPE dominance. Historically, 
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the success of OPEP motivated many attempts to derive a 


(two-pion exchange potential). However, the 
successes of these initial attempts were limited, and their 
derivations were not entirely ambiguity-free. 

In principle, to calculate the 27 exchange contri- 
bution, one needs all irreducible diagrams which have two 
pion lines joining the nucleons. Because of the infinite 
number of these, such a calculation is impossible. There- 
fore, in practice, one is forced to neglect certain diagrams. 
In fact, if one neglects radiative corrections, one is in 
fact neglecting what could possibly be sizeable contributions. 

Some of the many papers published on TPEP are those 


of Sugawara and Okuboe 's, Gammel et nile tee and Partovi 


and Temanncele soe 


When used with OPE models, these TPEP — 
represent some improvement over OPEP; however, there are 
conceptual difficulties involved which make such calcula- 
tions unreliable to a certain extent ‘2? , The basic diffi- 
culty, that of not being able to take into account all 
possible diagrams (or, equivalently, the problem with any 
perturbation theoretic approach) sets a limit on the accuracy 
of the models. In any case, TPEP generally suffers from 
problems in treating NN pairs and in treating the meson- 
meson interaction ‘1?). 

Chronologically, after the early TPEP attempts 


(16) (one-boson exchange) 


came the introduction of OBE 
models. These models assume that the NN force is meson- 


mediated, and that the exchanged systems may be adequately 
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represented by the meson resonances observed experimentally. 
This is based on the fact that in strong interactions, two 
or more particles as a group, tend to behave most of the 
time as a single particle with definite mass, width and 
intrinsic quantum numbers. The hope here, then, is that 
uncorrelated multipion contributions may be neglected, or 
at least represented by equivalent resonances. One advan- 
tage of OBEP models over TPEP models is that in the latter 
case, such contributions as 31 exchanges can never be 
included; another is that calculating single particle 
exchanges is a lot easier than calculating two-particle 
exchanges. 

The type of heavy boson exchanged is restricted by 
the symmetries of the interaction e.g. they must have zero 
strangeness. One difference between the various OBE models 
is in the particular meson exchanged. Needless to say, 
the full list of heavy bosons known has not been exhausted 
by the various theories; by the same token, it should be 
mentioned that not all the particles that have been used, 
have in fact been experimentally observed. In addition, 
calculations are done assuming the bosons to be stable 
particles; however, this is not true. In short, OBEP 
suffers from uncertainties as to which bosons are "really 
there", and with regards to treating massive composite 
systems (e.g. p resonance) which decay rapidly, as single 


particles ‘17), 
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There are, of course, certain nice results from 
this approach, e.g. the prediction of a large spin-orbit 
force, and of a repulsive core (see ref. 2 for further 
references). However, these are only qualitative results. 
Among the OBEP can be mentioned those of Bryan, et eh ee 


(19) (20), 


Green and collaborators , and the Japanese group 


An extensive list is tabulated in ref. 2, and a review 
of OBEP has also been done by Ogawa ‘21? | 

This, then, is a cursory glance at the situation; 
though brief, it shows that despite the great deal of work 
that has gone into the 2N problem, and despite the great 
strides that have been made, the problem is far from being 
solved. In this work, a fresh attempt is made to construct 
a field-theoretical model to describe the N-N system, free 
of divergences, and independent of perturbation theory. 
It will be seen that intrinsically non-local, and explicitly 
energy-dependent potentials are in fact necessary for a 
proper description of N-N phenomena. Although the equations 
involved are exceedingly complicated, approximations can 


be made to enable one to extract information from them. 
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III MATHEMATICAL FORMALISM 


In this work, the interaction between the two 
nucleons is assumed to be carried via the exchange of 


mesons (scalar and pseudoscalar, in the two cases to be 


dealt with); this amounts to prescribing a field-theoretical 


model. The model is then quantised a la Capri ‘22), In 
this section, the method of quantisation will be outlined. 
In essence, this will be a resumé of ref. 22. 

In general, a relativistic quantised field w will 
consist of both creation and annihilation operators, due 
to the presence of both mass hyperboloids (or positive 
and negative energies) associated with the field. It is 
not just the relativistic nature of the field that is 
responsible for this - if one has a non-relativistic many- 
body problem, the same situation arises, because of ex- 
citations above the Fermi level, which leave holes behind. 
On the other hand, in a non-relativistic few-body problem, 
this is not so, and one can take fi@inls to be a pure anni- 
hilation (creation) operator. It is precisely this feature 
of the non-relativistic problem that facilitates the 
quantisation procedure. In both the models to be dealt 
with, the nucleons are non-relativistic, and therefore 


the technique becomes directly applicable. 


Consider the general class of equations of the type 
A) y2 3 >> ach ees > 
Case ir) lee a) = | d>xtdt'V(t-t';x-x') v(t, x) W(t" x") v(t! 
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This is, of course, a Schrodinger-type eqn. for the field 
wv, which is both a non-local, and a non-linear equation. 
One first assumes that there exists a complete set of 
physical states (in the domain of vv), viz. the vacuum 
|O0>, one-particle states with definite momentum lk>, two- 
particle scattering states of definite momenta ee ae 
two-particle bound states labelled by the binding energy 
En and centre of mass momentum Pp, (ee p>, and so on. 
Furthermore, one assumes equal time commutation (or anti- 


commutation) relations, defined between pairs of physical 


states; for the vacuum and one-particle states, this yields 


<0] [y(t,%) ,W (t,y)1,/0> = 6 (%-9) (zr12) 


& (k 


<k, | tv (t,) 0" (t,9)1, 1k5> “k,) 6 (Ry) (III.3) 


i 
The general formalism can be extended in a straightforward 
manner to n-particle states. But as one is primarily 
interested, in this case at least, in 2-nucleon systems, 
it is sufficient to go only as far as n = 2. The objects 
of interest are not the field yp or the states [Reeekees 


but rather the matrix elements of products of the field 


operators between these physical states; these matrix 


elements are then interpreted as the physical wavefunctions. 
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III.1 One Particle States 


In order to obtain an equation for the two-particle 
wavefunction, it is necessary first to obtain the one- 
particle equation. The one-particle wavefunction is 


defined as follows: 
uy, (t,x) = <O0|y (t,x) |k> (arya) 


Quantisation amounts to taking the expectation value of 
eqn. (III.1) between an n-1 particle state, and an n- 
particle state, thereby yielding an equation for the n- 
particle wavefunction. In this case, n=l, and the one- 
particle equation is obtained by taking the expectation 
value of eqn. (III.1) between |0> and |p>. One then 
inserts complete sets of states between the field operators 
on the righthand side of (III.1). However, peace ww) 
is purely annihilation (creation), the only intermediate 
states that contribute will be one-particle states (between 


b (t,x) and i om eauly and the vacuum (between pe (ese) and 


W(t',x")). Consequently, the equation for uy, (t,X) becomes 
us 4 EE gee ee “s 
(ing + 5p) Uy, (tx) = dwqdex "dtr (t=-tlsx=x2 ) u, (t,x) x 


P ue" Ka (eK) (III.5) 


It is obvious from (III.5) that although the original 


field equation (III.1) was non-linear, because the field 


16 


| —— 
at ona tovew. oo Sea 


if 


es) z oe a a 


a 


h Us Pats ba a i = 
eee Cie he ee ee DP eee ee , 7 
Or $aID 27248 Og CML NES oF - dita iis OLS 

mt 


ste 


; — 
, 85678, efoisisg’ ree ts vagina” tt 
‘ ‘3 >t A ay ‘ A 


mo eid fos \t-t, 9865 Skat at “.aoLsae 
ae Piphta ya bermtbido @k 
“Gi B30] Aeswised ({ 37 
KA S& rE TR 5 G -4 Se 


. $Y gevenvsd ,rsvewoH  ~(2DV2EL) 26 ob be brsdatig t a 
cE y ine Sri : (nrotyses%tD) FI OT mea 4 rite pip 
i 


od) petete olotitasd—sno sad CPkw ssudkrdat ws 2 ile 2a3s 


a 
. 


s : Pe, 
& ff | nf at ‘ 7 : 4 : see | a ée) en a J Pais 
Pisit edt sex ai ,wsonrti-non esw (i .Trt) noLte pe bf 


he the tn 


= 


_ - : y 


; £ 7 


a 


foeneither only creation or only annihilation, the one- 
particle wavefunction does not couple to any other wave- 
BUNCE LON. 

Upon making the Ansatz 


d >> 
* -i (E, t-k-x) 
u, (t,x) = Ae ; CL LTE6)) 


eqn. (III.5) is converted into an integral equation for Eye 
which can be written symbolically as 


k2 
Ex = =u t (27) V, CIEL) 


it 2S, Of Course, not at all obvious from (IlI27) that 
this is an integral equation; in fact the energy just 
seems to have been shifted by an amount proportional to 


Vi° However, ve is actually given by 


~~ 


A 3 oi ne Pei a 
ve | d-g VIE, Bik q] (TII<.8) 
with 
Sy ; > > 
V(w,k) = —~, fat demay (sd ee ee 
(27) 


(Tier ¢9)) 


Therefore, as Vee itself is an integral involving Exe eqn. 


(III.7) is obviously a non-linear integral equation for Exe 
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(x * 1-3) t= 


3(xk\3)V +b 26] y , = (A, w) ¥ Se 
(we) ; 
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uy; 


(Note also that the factor 27 in (III.7) arises from the 


use of the ETCR (or ETACR), which gives ac = (dn) >.) 


Instead of actually solving for Ex from (ils. 7s; 
which is a formidable task, what one does is to insert 


the "free" value of the energy i.e. 


Pg 
k 2M 
into the expression for Vie Then Vi. is reduced to purely 
an integral term; as it is the difference between the 
"total" and "free" energies for a one-particle system, 
it is identified as the self-energy. It should be mentioned, 
at this juncture, that this definition of self-energy is 
consistent with the corresponding situation for the static 
case, where V(t-t';x-x') = S(t-t')V, (X-x"), which has also 
been discussed by Capri ‘22? 
To conclude, the one-particle wavefunctions are, 
in fact plane wave solutions, with the energy being defined 
by eqn. (III.7). Since the a Cf, x) form a complete set, 


one need not look for other solutions to the one-particle 


equation. 
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III.2 Two Particle States 


Because of the two possible types of states that 
° s ° . =A ai 
could occur in this case, viz. scattering states IP} Po” 
and bound states |E., D>, one naturally has two types of 
wavefunctions, and consequently two equations governing 
their behaviour. Scattering and bound state wavefunctions 


are defined as follows: 


<0) (t,2) (t,x) |pyB5> (eet 0) 


a 
a 
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One can derive, as in ref. 22, equations for the 
: >> >> 54; 
objects F p (t,xX,zZ) and Sap URIS However, the equations 
oa? 
one eventually wants are for functions of simpler structure. 


These functions are obtained as follows: as one is inter- 


ested in stationary states, one writes 
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From this point on, the one-particle energy E. is assumed 
to be just the free part, i.e. self-energy is neglected. 
Then, the energies Epi Po) and E,, (p) are written as 
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E, being the binding energy of the nth bound state. Further, 


the centre of mass motion may be separated out as follows 


= 
ot > X+ 
i(p,+p.,) ° (—5— 
. (xz) ‘=e 24 Cee (seo) (eG) 
FPO Pi RPS 


with a similar relation for the bound state. The equations 
one eventually obtains are for the functions £ (©) where q 
and r are related to the relative momentum and relative co- 
ordinate variables respectively. Substituting eqns. (III.12) 


to (III.16) into the equations for F and F__, one is led, 
Eno pn 


after fairly simple but tedious algebra, to the following 
pair of Schrodinger equations, in the centre of momentum | 


frame: 
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The V's are given by: 
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(The above equations differ from the corresponding ones in 
ref. (22), because the latter contains errors.) In addition, 
if one works out eqn. (III.3), it can be seen that the ETCR 


(or ETACR) reduces to the following condition: 


| a3g £ (F)f (e') + J£. (FE. 2") = —E Qs @-E") #6 (f42')] 
q q an n (27)? 
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What has been done here is the following: starting 
from the original equation for the field wp, equations have 
been derived for the possible two-particle states. Equations 
(id Ind gatandy(Iill.18) form a pair of coupled, non-linear 
Schrodinger equations which describes the wavefunctions for 
the stationary states of two particles interacting via 


V(t-t';x-x"). However, the coupling is only between the 
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two-particle states themselves, and between the two-particle 
and one-particle states through the term containing We (tat) 
i.e. there is no coupling to higher-particle wavefunctions. 
A certain amount of decoupling has therefore taken place. 
In addition, the terms in the square brackets in eqns. (III.17) 
and (III.18) are interpreted as "potentials" within the 
Schrodinger equation, and are, in general, explicitly non- 
local and energy-dependent. This is purely a consequence of 
the original equation of motion, and if one has a model for 
the N-N interaction governed by such an equation, one will be 
led quite naturally to non-local, energy-dependent potentials. 
The solutions eglsl and ale) are constrained further 
by requirement (III.24); it is a completeness condition on 
the two-particle solutions, which is, in this form, a gee. 
malisation condition. 
One point that has to be mentioned at this stage has 


to do with the self-energy term. One has identified the term 
27M | dive (per ete) 
in eqns. (III.17) and (III.18) with self-energy. This iden- 
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tification becomes obvious when one sees that V, (x) and Vi. 
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Before calculations can be done, the self-energy has to 
be explicitly subtracted out. This will be discussed 
further in the next section. 

What this technique does, therefore, is the fol- 
lowing: starting from a non-relativistic field equation 
of the type (III.1), one derives a set of coupled, non- 
linear n-particle Schrodinger equations. Because of the 
non-relativistic nature of the model, these equations can 
be partially decoupled, though they remain non-linear. 
When, therefore, one approaches the N-N problem, one will 
be led to equations describing the wavefunctions of two- 
nucleon states, in both scattering and bound state channels. 
From these equations, physical properties of the N-N system, 


such as phase shifts, binding energies, effective ranges, 
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etc. can, in principle, be extracted. As they stand, however, 


(III.17) and (III.18) are obviously intractable, and they will 


first have to be simplified. 

The two models to be considered are the scalar and 
the pseudoscalar models. In both these models, the equation 
for the nucleon field can be reduced to the form as given 
by (III.1), and so the technique described above goes 


through in a straightforward manner. 
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IV THE SCALAR MODEL 


Field theory describes the interaction between 
nucleons as being transmitted via particle exchange. 
Symmetry and invariance requirements dictate the particle 
exchanges possible. In the two models to be considered, 
it is assumed that meson fields act as the intermediary 
between the nucleons. The difference between the two models 
is in the type of meson exchanged. In the scalar model 2) 
there is no spin or charge; obviously, such a model does not 
give a true description of the physics. In fact, "truly" 
scalar mesons, which are described by equations derivable 
from a Lagrangian, give a repulsive potential. Therefore, 
in order to obtain an attractive potential, and hence simu- 
late pseudoscalar mesons, one uses, in fact, field equations 
not derivable from a Lagrangian. The pseudoscalar model, 
on the other hand, is more "physical", and does include spin 
and isospin. However, because of the far more complicated 
structure of the pseudoscalar theory, the scalar model is 
considered first. In this way, one is able to develop 
more easily, the techniques required, and also to gain 
insight into the computational problems that exist, and the 
complications that may arise. 


The equations for the coupled nucleon and meson fields 


in the scalar theory are 
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2 ee - > > 
(UO + m’)o(t,x) = - gp (t,x) p (t,x). (IV. 2) 


Here, ~, » are the nucleon and meson fields respectively, 
and g is the coupling constant. Eqn. (IV.2) can be solved, 
at least formally, and the solution written in integral 
form as follows: 

(t,x) =-g | Bleue ew Ce we x jdetd e, 


(IV.3) 


There is no inhomogeneous term in (IV.3). One reason for 
this is that one assumes that there are no free pions - 
all the pions created are absorbed. Furthermore, our 
interest is restricted mainly to elastic N-N interactions. 
In fact, this is also the reason for the choice of A as 
the appropriate Green's function. Eqn. (IV.3), when sub- 


stituted into (IV.1) gives an equation in wp alone; viz. 


2 >> > > 
(ise yaa) W(t) = - g° | A (t-t';x-%") (t,x) (t", x") * 


x w(t',x") dt! aex' . (IV. 4) 
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This has been achieved because the "source" term for the 
pion field is dependent only upon the nucleon field. 
Although (IV.4) is non-linear, it is independent (at least 


explicitly) of the meson field, and is of the type given 


by eqns (Iils1)..., The, potential, in this case, is ee 
> > 
i (wt-ke x) 
V(t,x) = - g7A(t,x) = i 218 
pene w FG x) = —4z |dw ak =——— 
(27) wo -k“ —m 

Therefore, 

V(w,k) = SiN el (IV. 6) 


In this model, physical states are defined, as outlined in 


Section III. Therefore, one-particle states are defined as 
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and this leads to the integral equation for EX: 
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The ETCR have been used to obtain lal? = (Qny7>s ie 


second term in (IV.8) is what one refers to as the self- 


energy. 


Iv.l The N-N Equations 


In the 2-nucleon system, there is only one bound 
state; namely, the deuteron. So, the label n used in 
Section III (for the number of bound states) will be 
dropped, as n= 1. The scattering states, therefore, 


> 


are of the type [K, oa and the bound state is |E k>. 
As in Section III, one can obtain equations for the cor- 
responding wavefunctions. They will be identically the 

same as eqns. (III.17) and (III.18), with the structure 
Omcne V's given by eqn. (IV.6). 

At this point, one makes an additional assumption - 
one assumes that the contribution to the non-local potentials 
from the snpermaad tee scattering states is dominant, and 
that the bound state contributions may be neglected. One 
reason for this is the following: in the pseudoscalar 
model, it will be seen that in the equation for the deuteron 
wavefunction, the only intermediate states that appear are 
scattering states, i.e. the bound state does not couple to 
itself. Therefore, as one is interested in the scalar model 
mainly as a guide to handling the more physical, pseudo- 


scalar model, one is quite justified in dropping the bound 
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state contribution in the former. There is also another 
reason: if one solves the Schrodinger equation with a 
6-function (or zero-width) potential, for both scattering 
states and bound states, one can then evaluate integrals 


of the type 
3 hee > —* ame 
| d~q Vag (t ue ee ae Ce ) 


where the Supe are these solutions. (Vee stands for the 
possible V's as tabulated in (III.19) COmsGid Le Vo eu) "WE 

can then be seen that contributions from intermediate 

bound states are small compared to those from the scattering 
states. 


Upon dropping these terms, one gets the following: 
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For the energy range one is interested in, where p 
is small, these potentials are almost identical i.e. for 
low energies, both types of states see the same potential. 
In addition, (IV.10) and (IV.11) are related through the 


substitution 


There dhe two points to be taken into consideration 
here - firstly, eqns. (IV.10) and (IV.11) are non-linear in 
the wavefunction, and therefore difficult to solve. The way 
this is achieved here is, in a sense, iterative; we only 
carry out this iteration to first order, and so the inter- 
mediate states are assumed to be free states. Therefore, 
properly normalised and symmetrised plane wave solutions 
are inserted into the integrals in (IV.10) and (IV.11). 
Hence, the equations have been linearised. The second point 
has to do with the self-energy term, which is still present 


in the equations. Once again, one uses the static model, 
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which has been solved by Capri, as the guideline. Twice 
the one-particle self energy term, i.e. a term 
47M | d>r' V(e-r') £(r") is subtracted from the right hand 
side of the equations. It will become apparent later, that 
Sucre a subtraction is necessary for the equations to be 
consistent. 

When these points are taken into consideration, the 
equations become 
5) 
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where the superscript stands for free solutions. These 
equations will be examined in both configuration-space and 


momentum-space. 
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ive 2 CONLTgGUuLacron Space Potentials 


The ground state of the deuteron is a spin triplet, 
isospin singlet state. Therefore, for the wavefunction to 
stay odd under particle interchange, the spatial part has 
to be even under parity. Therefore the intermediate wave- 
functions have to be symmetric in the spatial coordinate. 


So we write 


¥ 
Y 


cae ts) = CletP 4 e7tP't; 3 (IV.16) 


where C is determined from the ETCR (eqn. III.24), and 


turns out to be 
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Non-local potentials Vaeetiy o) and Ve" (f-£") are defined 


as follows: 
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The superscript NL stands for "non-local"; with these 


definitions, (IV.14) and (IV.15) become 


Oe ae) (IV.20) 
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Eqns. (IV.20) and (IV.21) are just Schrodinger equations 
with non-local, energy-dependent potentials which describe 
the scattering and bound states of the 2-nucleon system. 
It would, of course, be very desirable to have the Vibig 
explicitly as functions of r and r'; however, the actual 


task is formidable. Consider first, the object Q,(E,r") 


which is 
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and 
Vy (2,0") = MQ (r,r') - 2mM Vo (e,£'), 
Onna") being analogous to Q,(£,5"). 


Using the results of Appendix A, where these integrals are 


evaluated, one gets that 
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Because of the highly complicated structure of the arguments 
within the exponentials in eqn. (IV.23), the k-integral 
cannot be evaluated in closed form analytically - one can 
look at the problem equally well in momentum space, where 
equations for the T-matrix may be derived. As physical 
quantities are easily obtained from the T-matrix, the con- 


figuration space potentials will be left as in eqn. (IV.23). 
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IV.3 Momentum Space Considerations 


Eqn. (IV.14), the equation for the scattering wave- 
function fa) reads: 
3 


(Ve+p*) £, (F) = - 27M | d°r'! VE Geez”) 
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In momentum space, the corresponding wavefunctions are 


defined by: 
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The equations for Vy (K) and Vy (K) are obtained in a straight- 
forward manner from the configuration space equations 


(Appendix B). For V (kK), in fact, the equation is 
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To obtain the free intermediate wavefunctions 
gah: one uses eqns. (1IV.24), (IV.16) and (IV.17); 


these give 
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Substetutang (1Ve27) into (1V.26)° yields 
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the derivation of which is in App. B. One can now combine 
the coefficients of Vv, (K) accordingly, and with a function 


F(p,k) defined as 


quapti52662 HRESe 29 
F(p,k) = 27M f d>q Gj faPene 200, gal] -G[ pe get (IV.29) 
eqn. (IV.28) becomes: 
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[p?-k*-F (p,k) 1 (Kk) = 20m | d>q p/P ark y (q). 
Pp M Pp 
(IV. 30) 


This is, therefore the integral equation for the wavefunction 
v, (k) in momentum space (analogous to the integrodifferential 
eqn. (1V.20)) in configuration space). oF (p,k) containsythe 


self-energy, and one sees that when k=p, from eqn. (IV.29), 
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F(p,k) vanishes. The relevance of this becomes obvious 
in Section IV.5, where the self-energy subtraction is 
examined in closer detail. 

One object of transforming to momentum space is to 
enable one to obtain an expression for the T-matrix. In 


terms of the wavefunction vk) , the full T-matrix is 


defined by !?>) 
VO (kK) =) es 6 (DER) +6 (p-k b heee t(k,p,p°)+T(-k,p,p*) 
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(Iv. 31) 


The *T*s in eqn. (IV.31)' are half=off shell quantities, ‘the 
label pe indicating that p is on-shell. From now on, the 
p? will be suppressed, and we shall write the T-matrix as 
T(K,p). 

Experimental quantities, like phase shifts, etc., 
are related to the on-shell part of the T-matrix, viz. 
T(p,p) (2°). Substituting (IV.31) into (IV.30) gives the 


equation for the T-matrix: 
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with 


The equation for any given partial wave may be 
obtained from eqn. (IV.32); however, it is easier, in fact, 
to do this using the equation for V, (K) (eqn. (IV.30)) (see 


Appendix C). Eqn. (C16) for the s-wave T-matrix, To is: 
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can obtain a great deal of information from the on-shell 

Dale ot tea it is obvious from (IV.33) that one will never 
have the on-shell part until one has solved (IV.33) i.e. one 
needs the full Te before physical quantities can be extracted. 
However, by going into momentum space, one has reduced the 


integrodifferential eqn. (IV.10) into a purely integral one, 
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PV C45 {Born Approximation 


The most frequently encountered approximation for 


ee Eye 


solving the Schrodinger equation is that of Born 
amounts to a series expansion in "powers" of the interaction 
potential. The first term in such an expansion is referred 


to as the Born approximation. If one considers a Yukawa 


potential of the type 


Pl eel. L eh craw te ok! (IV. 34) 
(27)2m |m2+(Ak)2 m+ (zk)? : 


Ln WVe34)%, Ek and Ak are the total momentum and momentum 
transfer respectively, and the numerical factor comes from 
normalisation. 

In our case, the Born term can be derived in analogous 
fashion (i.e. from the inhomogeneous term in the equation for 


the full (non-partial wave decomposed) T-matrix (eqn. IV.32)). 


=> 
Then, one goes on to the mass shell by making || i 
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Dropping the self energy, this becomes 
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Comparing (IV.34) and (IV.36), we see that the Born approxi- 
mation for the Yukawa potential is reproduced with the 


coupling constant given by the relation 


V 2 


[Ome ge? ie 
Wom | ame 
2 V, 1 

oe — = — mg) e (IV. Sy) 
(27) 27 


For the Yukawa potential of the type we have considered, it 


is possible to obtain the deuteron properties (e.g. binding 


) (28) 
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Substituting these values into (IV.37) gives 


= O04 (IV. 38) 
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In our bound state calculation (to be discussed in 
Section IV.6), a variational approach was used to obtain 


a maximum value for the binding energy E Therefore the 


b° 
corresponding coupling constant obtained will in fact be 


an upper limit to its true value, and hence corresponds 


to stronger binding. The value we obtained was ae = 4.6, 
Or 
2 
—i_, = 0.0185 
(27) 


This value is in fact larger than that obtained from the 
Yukawa potential, viz. 0.0164, but at the same time is close 
enough to it to indicate good agreement between this model 
and the Yukawa model, at least within the Born approximation. 
At the same time, these two values seem to indicate internal 


consistency of the model. 


IV.o Sei ieinergy Subtraction 


In Section IV.1, the self-energy term was taken care 
Ofein) aavery particular way. §ln fact, cxactly stwice -tnesone— 
particle self energy term that appeared in the original 
equation, viz. (IV.10), was subtracted out. In this Section, 
it is shown that this is exactly what has to be done in 
fact, for the equation for the T-matrix to be free of incon- 
sistencies. In order to show this, one considers eqn. (C16) 


for the s-wave part of T, which is 
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If one assumes that a quantity which is a x the self-energy 
is to be subtracted i.e. if, after subtraction, eqn. (IV.10) 
becomes 
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This becomes obvious when one looks at (IV.28). So the 


quantity F(p,k) as defined by (IV.29) is modified to read 
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If now, one considers the on-shell limit of eqn. (IV.39), 
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which implies that, since the on-shell T-matrix can never 


be zero for all momenta, 
F(k,k) = 0 (IV. 42) 


Eqns. (IV.41) and (IV.42) imply that 
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So one is led to the conclusion that if V,(f-r") is to be 
interpreted as self-energy, then a = 47M, and in fact twice 
the one-particle self energy has to be subtracted. As 


mentioned previously, this is consistent with the subtraction 


in the static case, where it is obvious. It amounts to 
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subtracting twice the single-particle self-energy, as it 


is a two-particle system. 


IV.6 The Bound State 


In Section IV.1, it was pointed out that the bound 
state potential is obtained from the scattering state po- 


tential by the replacement 


The s-wave bound state equation can therefore be written 


down quite easily, by using eqn. (C9). One first defines 
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in analogy with (C1), and (C6) respectively. Then the re- 


quired equation becomes 


[k74ME+F, (k) 106°) (k) = (Iv.45) 


b 


Pie 


De RN ODL 8D ng ee 
(gq eu 4m“ (q-k) “-4M“m 1) a) 


2 ro) 
Mg 
7 | ood: ERY. NY: 
0 (q°+k°+2ME) “-4M" (qt+k) “-4M"m 


Meas 


42 


a? v 
' 7 tf 
i “e 
" 
i 
7 
S. rm 
Bru sod edz 
p Ls ’ are 
e= oe ey of 
Jy sate pilresse st 
‘@ 
nettitw sd eioleren: 
pas al 
Bamtieb se1tlt eno 
. 7 
! 
~ ym *\ 4 _ 
2 «V1) A+p 
' i. 
P h tr , D) F 
\ itd d Lt) { Al | } “sf v 
ae 
= <tr 
~9t ed neAT . 
' 
yy — 
ic®.Vi 
% 
A ; _— 
( ‘ (0), Mh ~ 
Pay ee 
x 7 > _- 
ua 


toads 3 se batntog esw +) 


‘ery moat = indo 


. (#9) .ape paheu yd .yitese os tup'g 


+r eudasaanes (180s ‘bas (197 sot ywotane 


R S (4t-p) * 2 ab-* Lt; cass ae® D Y las pbp. 
Mee” (A4p) "MB" ¢ ans + a+” Pri . 


aion | 
4, a aor 
A ity tem 


4is wn 
i 


Mota Hoomene cae 


pe 
~ : ‘ ‘ 
i oe 
nso nobdeupe otd29: buwod. ove 


a a) 


v4 
Rie : i 
g° 
D 
ist aise reat. 
(i, p) Vi shite = (0) 2) 
se 


_ 


| a 
: 7 iid Sate an 
; Wate a 


7 


. 
= «ay wrt, grrans 


Ld 
’ 


7) A ak a. us 
oo ee ; 
i a | 
is 


OL, 


(0) (y) 4 om Ky?) (x) - Fy, () ¥(°? (x) 


2 ve 
2 


eee Dane. 2 
(q*+k*+2ME) “-4M* (q-k)“-4M“m*|_, (0) 


» fore) 
0 (q +k°+2ME) “-4M" (qtk) “-4M"m 


2 (21) *k 


(IV. 46) 


In principle, one would use a variational technique 
to determine the binding energy E, given the coupling con- 
stant ace However, what is done here is to determine the 
value of i for which E is both a maximum, and also close 


to the physical binding energy. So one needs to maximise 
ME = - 4} yl) (Ky |2ax +] k2R_ (ck) [w8° (K) |2ak + 
0 b b b 


(q?+k*+2ME) 2-4m* (g-k)2-4m?m*| | 


2 roe) oe) 
st | kak| qdq an|-S.-*, 55 o5 
(q +k" +2ME) “-4M™ (gtk) “-4M"m 
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x 960) 1H) ar | + 
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One form of the deuteron wavefunction which has 
Simple analytic properties, and which closely approximates 


the square well solution is the so-called Hulthén wave- 


function ALAR 
f(r) = Cle t/%-e 7/8) (IV. 48) 
where the parameters are experimentally determined. (IV.48) 


is a difference of exponentials, and its fourier transform 
would be a sensible choice for BAe eK); the trial wave 
£uncrevonein (1V.47). “Chere is -alsimpler solution, viz. the 
solution to a square well potential which is a single ex- 


qe aes The increase in the degree of difficulty 


ponentia 
of computation of (IV.47) when one goes from a single ex- 
ponential to a function like (IV.48) is far too great to 
justify use of the Hulthén solution. 


e pL: . e e 
For an exponential of range 5 in configuration space, 


ve?) (k) is given by 


(IV.49) 


vn?) (kK) 
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With this, (IV.47) may be written as 


2 
- 2 kak Kk Fy, (kK) 
hea eS Bey yk ae ivr aes 
0 (k7+a7) 0 (k*+a*) 
ee aD DE in wey hae Gee 
, _Mg ; | ak f Aaa eid Sea ee aCe auans J 
0 0 (q7+k-+2ME)“—4M* (q+k) 7-4M7m 


k OS cl? OR 
Xx —— ; (TV .50) 
(ke40°) * eee | 


The three terms can be identified respectively as corres- 
ponding to kinetic energy (k), self-energy (s) and potential 


energy (p), and therefore (IV.50) is rewritten as 


with 
Dike ale aehe 
Tio= BAA | (Iv.51) 
Hs 0 (k2+a2) 4 


2 
o k FY (k) dk 


2 
Tosca | —5-- 
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(k2407)2 (q24a2)2 


Ty Can be evaluated exactly; in fact 


Tr. and Bt on the other hand, can be done analytically up 


to a stage; beyond that, the computation has to be numerically 


performed. 

The calculation performed is the following: for a 
given value of cn the expression (IV.50) is computed over 
a range of values of a; the maximum value of E, and its 
location in a are then obtained. This is repeated for 
different values of the coupling constant, until the best 
pair of values of a and g that gives an acceptable value 
of E is obtained, (0a was restricted to values below the 
one pion mass). It turns out that the computation is very 
sensitive to changes in the coupling constant, in the sense 
that the maximum changes value and location quite sharply 
with changes in g. The set of values that has been selected 

Z 


was at g” = 4.6, a = 122 MeV, for which the maximum is 


2.206 MeV. When g is increased, the maximum increases 
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and moves to a larger value of a. Whereas, for the square 

well, the range of the exponential is a (VME = 45.69 MeV), 
Vv ME 

the range one gets here is more of the order of one pion 

mass. Tables I - III contain the results of the calcula- 


tion for three values of ae Appendix E outlines the pro- 


cedure of integrating I, and Pia 


Iv.7 Evaluation of Phase Shifts 


The phase shifts for the &th partial wave are given 

by the corresponding on-shell T-matrix via the relation 
ge 
T, (p) = -o TDM e sind, (p) (IV.54) 

Here, aS in subsequent sections, the on-shell T-matrix will 
be labelled by the on-shell momentum i.e. T(p,p) = T(p). 
To obtain the phase shifts, therefore, one has to solve the 
integral equation (IV.33), and then go on to the mass-shell. 
Lt is quite obvious that an)analytic solution to (1V-.33), 
though desirable, is certainly not easily obtained. Therefore, 
the solution is obtained numerically using a technique of 
matrix inversion. 


For s-waves, we first rewrite (IV.33) as follows: 
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Table I 


Binding energy calculation (for ar = 4,25) 
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Table II 


Binding energy calculation (for aa = 4.60) 
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Table III 


Binding energy calculation (for g* = 5.00). 


0-550 Li Salar 

80 0.2024 2/0383 
90 0.2546 =O. L628 2.5497 
100 Oro LZ =O LIGU 3.0218 
110 Omsn Ly 002383 3.4247 
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with 


2 
A(k,p) = —_~ 2 M(p,k)V_(k,p) (IV.56) 
(27) P 
2 
B(k,p) = —_, 2 M(p,k) M (IV.57) 
(27) 
ee (q74+k?-2p7) 2-4m? (q-k) 2-4m2m? 
MN USISH Ness oesbe) Mere poy amma T Tare Dy (NY Sh) 
q (q°+k“-2p*) “-4M* (qt+k) “-4M“m 


Equation (IV.55) therefore becomes 


se q?v (k,q) 
To (kyp) + B(k,p) | dq +t, To (arp) = A(k,p) 
0 p -q +in 


(ase )) 


which is then written in the form of a matrix equation. 

p is kept as a parameter in the calculation, and for a 
given value of k, T (K+) is obtained for various values 

of k. In other words, half-off-shell matrix elements are 
obtained by this technique. The details of the computation 
are contained in Appendix D. 

The real and imaginary parts of T, (kp) are printed 
out, and the corresponding phase shift 6, (Pp) (both real and 
imaginary parts) as well, utilising eqn. (IV.54). In all 
the computations done, the imaginary part of the phase shift 


is at most of the order of 10713 radLansthelThus.1sato be 
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expected, of course, as one is dealing with effects 

below the one-pion threshold i.e. no absorptive or dis- 
Sipative effects enter. Also, for the same reason, uni- 
tarity can also be shown to hold almost exactly: the 
potential is real, and no approximations made should affect 
unitarity. In fact, the half-off-shell ae ie, condition 


reads 


Im T, (k,p) = - “F Dav enora dp) 

and the data show that for any value of p. and k, this is 
almost exactly valid. The difference between the two sides 
of the relation is of the order of yee this arises from 
the numerical (as opposed to the exact and analytical) 
approach, and is, for practical purposes, equal to zero. 

Three factors have to be considered in these calcu- 
lations. The computations can be done for various values 
of the coupling constant; in addition, sensitivity to the 
cut-off and to the self-energy can also be tested. The 
cut-off, if present, restricts the range of the q-inte- 
gration; ideally, one would like to see no variation of 
the results.with a.variation in. the,cut-off... The, larger 
the cut-off, one would expect more structure introduced 
into the integral in eqn. (1V.59); hence if the structure 
adversely affects the numerics, lowering the cut-off should 


reduce the structure, thereby improving the numerical work. 
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As for the self-energy, the matrix inversion can be per- 
formed both in the presence and the absence of the self- 
energy term. 

The most interesting results are shown in Graph I. 


This contains 3 curves, all with the coupling constant 


2 
Ss (OPO LSS 


(20) 7 


The three curves cover three possible situations with a 


fixed coupling constant, viz. 


a) Excluding self-energy; cut-off = 1 nucleon mass; 
b) Excluding self-energy; removing cut-off; 


c) Including self-energy; removing cut-off. 


It is obvious that in going from a) to c), more 
effects are being incorporated into the theory. As the 
curves indicate, the agreement between these three cases 
is reasonably good. The phase shift seems to drop from a 
value of 180°, therefore indicating the existence of a 
bound state. 

Further discussion of the results appears in Chapter 
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SCALAR MODEL: PHASE SHIFTS 
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V. THE PSEUDOSCALAR MODEL 


In reality, the proton and the neutron have dif- 
ferent charges - they form an isospin doublet; in addition, 
they are spin-% particles. As opposed to the scalar theory, 
these features are incorporated into the pseudoscalar 
model. Consequently the states can no longer be simply 
defined by momenta and binding energies; spins and iso- 
Spins have to be included, too. With all the indices sup- 


pressed for the moment, the equations of motion are: 


2 


cea tt on) v(x 57%) 


ox 
fe) 


i £ te oy (x, x) 898 (x, x) (v.1) 


(D+ my (x0,%) = - £5970" (x, ,X) 0 TH (x) (v.2) 
Here, VW, $ are the nucleon and pion fields respec- 

tively, 7 and o are the usual isospin and spin matrices, 

and fo is the pseudoscalar coupling constant. (The time 

variable is Xo, as opposed to t, which has been used until 

now, and which will now be reserved for isospin.) Eqns. 

(V.1) and (V.2) ‘are analogous to (iV.1) and (IV.2), and, 

as before, because the interaction term in (V.2) is $- 

independent, it can be solved for Oia) and inserted 

into (V.1) to give the following equation, in all its glory, 


for the nucleon field: 
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(As before, we have no free pion field, and hence only the 
principal value Green's function.) In (V.3), primed 
(unprimed) suffixes refer to isospin (spin) indices; there 
is summation over repeated indices. Apart from the compli- 
Cations, Of ‘spinwand? isospin,» (V2.3) editfers from (1V.4) in 
that it contains the double derivative 9+“. the field 
operator ~ is, of course, a matrix in spin-isospin space, 
as evidenced by the suffixes. 

As already mentioned, the physical states of this 
theory are no longer defined as before. The vacuum |0> 
stays unchanged; one particle states are defined by momentum, 
spin and isospin, and their third components. Therefore, 
in full, one should write such states as |pSTs,t,>, S(T) 
being the spin (isospin), and S3 and t. the corresponding 
third components. However, as S = T = ¥ in this case, they 
will therefore be suppressed; instead, the third components 
only will be specified, and will be relabelled s,t. So 
the one-particle states are |pst>. Two particle scattering 
states are, in analogous fashion, [pyPoSTs,t47 7 where now 
S(T) labels total spin (isospin), and S3 and t, are their 
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one bound state (the deuteron), it is given by | PE, STs, 3 


For discussion of the NN problem, this is as far as we 
need to go. 
Analogous to (101.2): and=(1l1L.3) , the commutators 


read: 
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With eqns. (V.3), (V.4) and (V.5), one has all the necessary 


equations. 


V.l1 One Nucleon Equation 


For the one-particle wavefunction, one tries the 


following Ansatz: 


<O|) 41 (Xo) |pst> = 6 66 


as x! tp (Xo *) (Vv. 6) 


UW, (XG r¥) should, ian fact, ‘be labelled by the total spin’ and 
isospin, which in this case, as mentioned before, can only 
be 4%, and are therefore suppressed. (The spatial part of 
the wavefunction is labelled only by total spin and iso- 


spin.) Once again, one assumes free solutions (which form 
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Application of (V.4) to (V.6) gives the value of A; i.e. 
one gets that 
oe fare 1 “iE xX ipex 
<O| Vi yt (Xr) |pst> i eka ore eu ate ce (v.8) 
(27) 
The equation for the one particle wavefunction is derived 


as before (App. F) and gives the following equation for 


> 
ha agree 
aie 
(inz - SOM a tot) in 


Z ' 3 ' 3 ' xs om pa tee en ‘ a Won 
3£6 | dx (oS GINS ho: Be hls reg Ens XA FX-X PL Tess posts ) 


ye) 


As was the case in the scalar model, this non-linear equation 
for u, leads to an integral equation for the energy, i.e. 
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In addition to these, one defines, for future use, a 
quantity We) ict. V(x) (TIT 719)}.) vas.-follows 
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W_(x) = | d-pe W (Vi 3) 
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So, as before, the true energy is the solution to 
(V.10). However, in practice, a closed solution cannot be 
obtained. As for the previous model, therefore, the self- 
energy will be dropped wherever E. appears in the kernel 


for the two-particle equation. 


V.2 The Two-Nucleon System 


A typical 2-particle wavefunction looks like 
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The equation governing the behaviour of objects such as 
(V.14) is derived exactly as outlined in Section III, 
once the spin-isospin dependence is properly accounted 


for. This is done as follows: 
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which serves to define the wavefunction of the two-particle 


ST > oe ‘ ; : 
F (x_,X,Z). A similar relationship exists 
tle 


for the bound state. One can then derive equations for 


fe (32 XZ) and For (x /X,2) as before (see App. G). Eqn. 
PyP>o 0 Ep “Oo 


(G8) is the scattering state wavefunction when S(T) is the 


system, viz. 


total spin (isospin). Because of the different products of 
spin and isospin wavefunctions (e.g. <s}v|Ss,><thv'|Tt,>), 
it is not apparent whether or not S' and T' (the indices 

for the intermediate states) can take all values. Rewriting 


(G8) in full gives: 
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It is obvious from (V.16) that each state labelled 
by (S,T) couples to a number of possible states (S',T'). 
It will be only by looking at each pair (S,T) that one can 
determine which values (S',T') actually do appear. Because 
the two particle states have both spin and isospin equal 


to 4, S and T can only take the values 0 and 1. There is 
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fairly lengthy algebra involved in actually obtaining 
the separate equations for each (S,T) pair, and this 
removal of the Clebsch-Gordan coefficients is done in 
Appendices H and I respectively. 

Since the deuteron bound state is an S = 1, T= 0 
state, and in fact, the only bound state of the system, 


certain intermediate state contributions to (V.16) will 
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drop out, i.e. the only bound state wavefunction is Fie 


and this will simplify some of the equations. 
Taking account of this, and using expressions (H9), 
(H10), (16) and (17), the following equations are obtained 


for the four scattering channels and the single bound state. 
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Eqns. (V.17) to (V.21) form the set of coupled, non- 
linear, and non-local equations which describes the 2- 
nucleon system in the framework of this theory. There is 
really no hope of being able to decouple these equations; 
therefore, they will be linearised as was done in the scalar 
model. It is also apparent that the non-local "potentials" 
which appear in these equations differ in all the equations - 
hence each state has an inherently different potential, 

each of which is both non-local and energy-dependent. It 
can also be seen that the presence of (5, °¥) (6,°¥) terms 
leads to tensor terms in some of the channels; in particular, 
in the bound state equation, where in fact, such a term is 
expected to be present. Also in the bound state equation, 
viz. (V.21) the intermediate states do not include the bound 
state itself. (This has already been mentioned in Section 
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As described in Section III (viz. eqns. (III.12) 
to (III.16)), these functions aie can be replaced by the 
centre-of-mass wavefunctions foe and equations (V.17) 
to (V.21), which are actually two-particle Schrodinger 
equations, can be reduced to single particle equations in 
the relative ‘coordinate. The actual derivation will not 
be done here, as the identical computation is to be found 


in ref. 22. In brief, however, it involves the following: 


first, consider the derivatives that appear within the 


kernel. 
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The equations that will be examined in closer detail 
are the scattering equation for the (0,1) channel, (V.17), 
and the deuteron equation, (V.21). Using eqns. (V.22) and 


(V.23), these may be rewritten as follows: 
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The corresponding equations for the other channels 
can also be brought to this form. Eventually, after de- 
fining the centre-of-mass wavefunctions (in analogy with 
eqns, (lit. 12) to, (111.16)))(. and after avcertain amount of 
algebra (identical to that done in ref. 22), one is led to 


the following set of equations: 
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In eqns. (V.27) to (V.31), one has the non-local 
Schrodinger equations governing the bound state and the 
various scattering channels. These are the pseudoscalar 
analogues of (III.17) and (III.18). Similarly, the poten- 
tials defined in (V.32) to (V.34) are analogous to (III.20), 


and (III.22). Once again, therefore, one has obtained 
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differential equations for the various states, starting 
from a pair of coupled field equations, (V.1) and (V.2). 
Most of our interest will focus on the scattering 
equation (V.27). It will be transformed into momentum 
Space, and, as with the scalar model, a half-off-shell 


T-matrix will be constructed. 


V.3 ' The S = °0,/..T ‘= 1, Channel 


Apart from the fact that this channel is simpler to 
handle, due to the absence of a tensor force, a great deal 
of information about neutron-proton scattering is contained 
in it. Therefore, we shall go through the calculations in 
some detail. 

When the self-energy is subtracted out (as in the 
scalar model), the equation for scattering in this Skee. 


viz. (V.27), becomes 
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where, in addition, the bound state contribution has been 
dropped for the reasons already discussed in the scalar 


model. The steps involved in going from this, the config- 
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uration space equation, to the corresponding momentum 
space equation, are exactly identical to those involved 
in going from (IV.14) to (IV.26) in the scalar model. 


This gives 
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In this case, however, there are two intermediate 
Stavteceto be considered, avi. the: (0)\imand ath) L)yestates, 
As the former is a spin triplet, isospin singlet state, 
whereas the latter is both a spin and isospin triplet, the 
spatial parts of the wavefunctions differ. In the (1,0) 
case, it will be even, while in the (1,1) state it will be 
odd. And since these intermediate state wavefunctions are 


taken to be free, one has (see App. J) 
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Eqn .u(Vii39) is identical in form™torAty.28)> analogous to 


the F(p,k) of scalar theory, one can define 
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It is obvious from (V.40) that on the mass shell, 
b= ky Glp,k). vanishes... This,s has, to, be so,,if,the equation 
(v.41) is to be valid (see Section IV.5). Furthermore, a 
T-matrix can be constructed as before (eqn. IV.31), as the 


symmetry is the same as it was in the scalar model (S = 0, 
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T = 1 implies the spatial part of the wavefunction is even). 


Therefore, 
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Substituting (V.42) into (V.41) gives the following equation 


for the T-matrix: 
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Once again the third label, pS, has been suppressed in T; 
the complete analogy between (V.43) and the corresponding 
scalar model equation, (IV.32), is very obvious. 

The equation for the s-wave part of the T-matrix 
may be obtained by partial-wave decomposing vp (k) i.e. 


eqn. (V.42). This gives: 
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The notation has been modified slightly in (V.44), in that 
the subscript o now stands for 2 = 0. The rest of the 


partial wave decomposition goes through as before (App. C), 
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which, upon defining N(p,k) analogously to M(p,k) i.e. 
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For this particular channel, therefore, eqn. (V.47) 


is the s-wave equation, and can be compared to egn. (IV.33). 


Although the integrals involved are more complicated, due 
to the structure of UL ey eths | the principle of solution is 
identical to that in the scalar model. There is no advan- 
tage to be had in discussing this in detail, except to 
mention that the approach is identical to that mentioned in 
App. D. 

For completeness, the equation for the T-matrix in 
the other three channels can also be derived, and are 


included here: 
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The corresponding partial-wave decomposed equations 
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Vote Born Approximation 


As in the scalar case, one can obtain the Born term 
in the S = 0, T = 1 channel. Once again, the self-energy 
is dropped, i.e. N(p,k) is made equal to l. Then, as 


before (eqn. (IV.35)), the Born term is 
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Eqn. (V.53) indicates that the Born term contains, in 
addition to an attractive Yukawa-type term, a constant term. 
In configuration space, this constant arises from a delta- 
function at the origin which, for local potentials, has no 
effect on the wavefunction. However, in our non-local case, 
the 6-function appears to get spread out, and since the 
sign of the term is positive, a repulsive core seems to be 


present. So the Born term appears to arise from a repulsive 
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core and an attractive Yukawa potential. This feature is 
quite different from the earlier model, which was purely 
attractive. One would expect the repulsion to "pull down" 


the phase shifts accordingly. 


V.5 The Deuteron 


The deuteron is predominantly in an s-state; however, 
the existence of a non-zero quadrupole moment, and the value 
of the magnetic moment indicate the presence of a d-wave 
admixture. This actually requires the presence of a tensor 
term in any potential which is supposed to describe the 
deuteron. 

Assuming that the deuteron wavefunction is a mixture 
of s- and d-waves, one can derive coupled equations for 
the two components. This leads to the Rarita-Schwinger 
equations '?9) first obtained in 1941. A similar approach 
can be taken in our case. The main difference will, of 
course, be that the potentials (central and tensor) will 
be non-local, and in principle, determined, whereas the 


potentials are usually inputs in the conventional treatment 


The deuteron wavefunction satisfies eqn. (V.31) which, 


upon subtracting out the self-energy, becomes 
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The term involving Wee will introduce a tensor term 


into the potential. (V.54) is first rewritten as 
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from the definitions of W and W In (V.56), the tensor 


operator is introduced via the relationship 
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So the (5, °K) (65°) term is replaced by a tensor term, and 
the (5, °5,) term adds to the central part of (V.56). 


Eventually, therefore, (V.55) becomes: 
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which may be written symbolically as 
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Eqn. (V.57) is a Schrodinger equation with non-local 
> > >> . 
central and tensor forces Volr,r') and Vip (xe ,x') respectively, 


which are defined as follows: 
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Eqn. (V.59) is now decomposed into partial waves. Since the 
deuteron is in a spin triplet state, (51°55) = 1. Following 
the method and notation of pepénedettit>*?, (where the 
analysis for purely local potentials is done), we can write 


the deuteron wavefunction as a combination of s- and d-wave 


components, i.e. 


as. 


where £,(r) and fy (r) are the s- and d-wave parts, and 
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usual spherical harmonics. If we assume that the deuteron 
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is in fact described by (V.61), we can decompose eqn. 


(V.57). In order to do this, we need the following: 
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Then we are led to the following coupled equations: 
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0 
Eqns. (V.63) and (V.64) are identical to the Rarita- 
Schwinger equations, except that the potentials are non- 
local (and inherently so), in this case. Both Vo (rer) 
and Vip (rer") may be written down, by just partial wave de- 
iets >> >> : 
composing the original Vo (rrr') and Vp(rer'), i.e. eqns. 
(v.58) and (V.59). So if,.as in the scattering channels, 


one uses free solutions for the intermediate states, one has, 
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in principle, a non-local potential describing the deuteron. 
Once again, however, the integrals cannot be evaluated in 
closed form. The point to note is that instead of stipu- 
Rating that Vo and Vin (in the local case) are, for example, 
square wells, or other simple potentials, we have, in fact, 
a potential derived from the theory. Examination of (V.58), 
for example, will show that, even in the local limit, there 
is a great deal of structure in the potentials. Although 
square well potentials are used for simplicity, it would 


seem that one actually needs potentials of a far more com- 


plicated nature to describe the deuteron. 
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PSEUDOSCALAR MODEL: SOME RESULTS 
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VI. DISCUSSION OF RESULTS AND CONCLUSIONS 


The technique that has been employed here had its 
roots in the solution of the anharmonic oscillator pro- 


plem ‘37? | Further investigations, particularly into the 


(22) indicated that this 


non-relativistic few-body problem 
was perhaps applicable to the interaction between two 
nucleons. Field theory has, of course, been used to tackle 
this problem before, and the vast quantity of literature 
attests to this fact. Most of the previous work, however, 
has been perturbation-theoretically inclined. 

The applicability of this particular scheme to the 
current problem depends upon two facts. Firstly, the 
equations of motion (or the Lagrangian, for that matter) 
are of a very particular type; the "source" term for one 
of the fields (the meson field in this case) is entirely 
independent of the field itself; i.e. there are no direct 
meson-meson interactions. This enables the fields to de- 
couple (at least formally). Secondly, the models are non- 
relativistic, and our interest is confined to (at most) 
two-particles. This few-body aspect enables the different 
particle-number states to partially decouple; i.e. for a 
solution of the two-body system, only one- and two-particle 
states need to be considered. 

We examined first, mainly for purposes of developing 
the methods required, a model in which the non-relativistic 


nucleons interact via the exchange of spinless, chargeless 
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mesons. This model, the scalar model, is far more 
tractable than a more physical one, not only because 
of the absence of spin and isospin, but also due to 
the simpler structure of the kernels that appear in 
the final Lippmann-Schwinger equations. 

By considering single particles in this theory, 
we saw that "free" wavefunctions, which form a complete 
set, do describe them. This implies, however, that the 


usual non-relativistic energy-momentum relationship 


can no longer be true. There is an additional term, an 
explicit form of which can only be obtained, in principle, 
on solving a non-linear integral equation for the energy. 
This extra term is interpreted as a self-energy; such an 
interpretation seems quite natural, at least for single 
particles, as there can be no other interactions present. 
Due to the complexity of the integral equation for Jaye we 
have not solved it explicitly; instead, we have identified 
the integral cao i.e. the difference between ES and the 
"free" energy a as the self-energy. 

In the two-particle case, equations are obtained 
with "potentials" which are both non-local and energy- 
dependent. There are only two channels possible in the 


scalar model, viz. scattering states and bound states. 


The scattering states were examined in both configuration 
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space and momentum space. In configuration space, an 
integral representation of the non-local potential was 
obtained; it is not possible to write down a closed, 

exact form of the potential in configuration space. From 
the momentum space equation, an equation for the half-off- 
shell T-matrix, T(k,p,p*); was derived. This is, again, 
an integral equation, and from the inhomogeneous term, the 
Born approximation was obtained. This term is identical, 
to within factors involving the coupling constant, to that 
obtained from a local Yukawa potential; in this approxi- 
mation, at least, the potential is attractive, due to our 
"strong-arming" the coupling constant. On shell, it is 

in fact a Yukawa potential. 

The s-wave T-matrix equation was examined further. 
Half-off-shell matrix elements were obtained; a lengthy 
algebraic and numerical procedure was required to do this. 
From the corresponding on-shell values, phase shifts were 
calculated. These phase shifts are plotted on Figure I, 
for three different combinations of parameters. The phase 
shift starts at 180°, and so indicates the presence of a 
bound state (the deuteron). After an initial rapid drop, 
the phase shift seems to level off at rather a high value 
(~ 60°). This would seem to indicate that there is not 
quite enough repulsion in this model. If our treatment 
of this model is to fit experimental data reasonably well, 


and if the phase shifts are to drop quicker (and, in fact, 
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take on negative values), repulsion would have to be 
incorporated in the form of a core. 

The results we obtain seem to indicate a good 
solution of the Lippmann-Schwinger equation; there is 
numerical consistency, and very little sensitivity to 
the value (or even the presence) of a cut-off. Figure 
I also shows that the inclusion of self-energy does not 
alter the results by much; this supports the contention 
that self-energy can, for all practical purposes, be 
neglected. 

The numerical results also satisfy unitarity very 
well. Any deviations from exact unitarity are far too 
small to be anything but numerical (roundoff) in origin. 
In any case, we would expect unitarity to Weticls as the 
potentials are all real. Also, the energies we are con- 
cerned with are low enough that only elastic processes 
are involved; therefore, the phase shifts should all be 
real. Numerically, the calculation yields values of the 
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radians for the imaginary part of 
the phase shift; this can be compared to values of the 
order. of..1 or 2 radians for the real part. 

As for the bound state in the scalar model, the 
treatment is slightly different. Using a trial wavefunction 
(a single exponential in this case), we look for a minimum 
in binding energy. Any possible refinement on the accuracy 


of the results by the use of a more complicated wavefunction 


(e.g. the Hulthen wavefunction) was far outweighed by the 
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added complexity of the calculation on doing so. The 
calculations exhibit a large degree of sensitivity to 
the value of que the coupling constant. By constraining 
the value of the parameter of the wavefunction (or 
equivalently, the range of the exponential) to be close 
tO vME, a value of ae = 4.6 was obtained. This was the 
value of g? used in the phase shift calculations, and 
also in the discussion of the Born term. 

One further point should be made regarding the 
scalar theory. The self-energy term was examined in 
both the zero and infinite momentum limits. In the former 
case, it showed that a renormalization of a few percent 
was all that was needed to correct for self-energy; on the 
other hand, the infinite renormalization is identically | 
zero. These calculations nave not been included; instead, 
the corresponding pseudoscalar calculations appear in 
Appendix F. 

The second model considered is the pseudoscalar one, 
where both spin and isospin are included. Apart from the 
additional matrix algebra required, the one-particle system 
is treated identically to that in the scalar model. Once 
again, a self-energy term can be derived. 

The two-nucleon system is quite a bit more compli- 
cated. There are four possible channels in which scattering 
takes place, and the single (deuteron) bound state. 


Equations here obtained for all five wavefunctions; also, 
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the T-matrix equations were obtained for the scattering 
channels. The potentials are more involved now, and in 
some channels (particularly the deuteron one), tensor 
forces appear, as they should. 
One of the simpler scattering channels is the 
S = 0, T= 1 channel. There is no tensor force present, 
and the method employed for the scalar model should carry 
through. But if one compares the equation for the T-matrix 
in this channel with the one in scalar theory, it can be 
seen that the kernel differs from the scalar one by a 
polynomial factor; also there are two (separately) divergent 
terms in the (0,1) channel which, together, will cancel out. 
Numerically, the results we obtained in this case 
are not as good as in the earlier case. Whereas, for the 
scalar case, numerical consistency was obtained using 
30 x 30 and 40 x 40 matrices, we had to use 60 x 60 matrices 
before a solution begins to appear (and sensitivity to the 
cutoff begins to vanish). This could arise from a number 
of things, one being increased structure in the T-matrix. 
In solving the equation, we have used linear interpolation 
of the T's; perhaps this is not enough to account for the 
structure. Taking a larger grid size circumvents this 
Gite rem ty. 
Using the physical coupling constant (£27 (20)? = 0.189) 
we did obtain a solution (see Fig. II); however, it still 


seems to contain too much structure. Also, it drops too 
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rapidly through zero. As opposed to scalar theory, therefore, 
there is obviously too much repulsion here. The calculations 
were repeated with various values of coupling constant, 
but this does not alter the degree of repulsion present. 
These difficulties seem to be tied in with a more 
general problem. Theoreticaily, it is exceedingly difficult 
to describe the s-wave interaction well. At short distances, 
exceedingly complicated particle exchanges come into play; 
this does not occur for higher partial waves. In fact, the 
problem can be seen quite clearly when one considers that 
the origin of the repulsive core came from the fact that 
the sx phase shift goes negative - a fact that cannot be 
explained otherwise. It is apparent, therefore, that the. 
difficulty of handling the 25% state adequately has been 
around for a while. In particular, the amount of repulsion 
needed is one aspect that has yet to be treated completely, 
georatiicaliyn> 
The two models therefore seem to differ basically 
in the amount of repulsion present. The Born terms on- 
shell, for example, differ basically by a constant. This 
constant, which is analogous, in the local case, to a 6- 
function potential at the origin,in this non-local case 
creates a sort of core, thereby causing the repulsion. 
With this in mind, a calculation was performed in the 
scalar theory, to see whether the ash data can be reproduced. 
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drops from about aoly However, once again, 84 stays non- 


negative up to at least 280 MeV. In addition to the phase 
shifts, the off-shell behaviour of T was examined. In 
Fig. IV, the matrix elements T, (p-k) (normalised to 1) 
are plotted for two values of EO? viz. 60 MeV and 180 MeV. 
Typically, these curves should go through zero; but 
because the phase shifts stay positive, so do these ratios. 
In the pseudoscalar model, this difficulty seems 
to be present only in the s-wave. This is primarily due 
to the constant term in the Born approximation mentioned 
in the previous paragraph. However, upon partial wave 
decomposing, this constant only contributes to the s-wave. 
Therefore, for higher partial waves, there will probably 
be improvement. A preliminary calculation, using only the 
Born term, seems to confirm this. 
At this stage, it is safe to say that we have a 
numerical technique that enables us to obtain a half-off- 


shell T-matrix, and corresponding phase shifts, froma 


theory such as this. However, only 2=0 has been considered, 


and though generalisation to g=2, etc. is straightforward, 
it is far from trivial. The structure of the kernels 
become far more complicated. Although the numerics appear 
to be acceptable, the physics in the model does not give 
very good agreement with experiments, at least for s-waves. 
But, as mentioned before, this probably stems from the 


nature of the s-wave interaction itself. Theory appears 
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to keep having a very difficult time with theory. 
Inherently, this technique is of value- qualitative 
features of the N-N interaction seem to arise naturally: 
a repulsive core; the tensor force in the deuteron; 
different potentials in the various channels, etc. 
Quantitatively, it enables one, in principle, to evaluate 
half-off-shell T-matrices. This in itself probably 
justifies further work, as lately, such calculations have 
been increasing considerably in number. 
As for the non-local potentials themselves (e.g. in 


the deuteron), both non-locality and energy-dependence 


appear naturally. However, the structure of these potentials 


lead us to believe that no local approximation is going to 
take account enough features of the original potential. 
Such approximations may be useful in very weakly or very 
strongly non-local potentials; here, the potentials are 
far too complex to be so categorised. 

For scattering, this technique should be applied to 
higher partial waves. If, as we speculate, the problems 
we have encountered arise from the nature of thes-wave 
itself, a calculation for %=2 would certainly produce better 
results. Also, it would be worthwhile to develop a method 
for solving all the four coupled equations. In principle, 
a (much) larger grid would permit this to be done. As for 


the bound state, the non-local potentials ie and Vin should 
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be examined more closely. The possibility of a numerical 
solution (in an iterative manner, perhaps) of the deuteron 
should be looked at. Whether one can obtain a reasonably 
good approximation to the (now) complicated potentials is 
a question well worth pursuing. There seems to be, at this 
stage, anyway, more work needed to see cerner (ears theory 


Willi stand or fall. 


APPENDIX A: EVALUATION OF THE POTENTIAL 


3 Be > > 
| dq V Bes 74° ¥ 


LS 
p 
> 
Seige dean. PP OLteD 
= (21) 4 f dq ak etrd Ye 2 v(-t pk +k] 
2 ah 83 ik. a*y 1 
[4s -4p“+k ene 
aM as a 
k 2 eae 
-i-° —_= Qq*y 
= 16mg” | ee” | Be ene a ips Se eS 
(4q -4p +k”) -16M"k°-16M'm 
(Al) 
Consider 
Di tqe=4 bank ic ae oMck 2 wee me 
2 7 * 44 73 
= 0 when Gay = NDS Sai eM noite es (A2) 


The various possible cases have to be considered separately, 


to determine the possible signs of q? 


Case (i) Ka Ze 


In this case, De ae 


oe, 


, 
M% Mi iL 
V 
S : " c 
(S48) “ee NM t q = ~p ° ers 
- 
| é ” — 
ylejsusgee betoebteamoo ed og eved aeeass sldieeog 
: y | ms 
c 
. .p 20 ampie eldtezog 
ia a 
t 
; a coy / 7 
¢ <A. 
i ee © Be 
lt ae 
a i we) is 


ay : ayswin 0 


i ' oh i ‘ : 7 | : ; 7 ra 
y ye 7 - _ - *& =f 


Se eee eS 


a 


But q* could change sign; 
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~ m2 (k24+m7) = 0 
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The energies one is interested in are below the one-pion 


threshold, i.e. p- < mM; in this range 


2 
gM? + 5p) < om /l+ P+ SB, 
M 4M 
so the only positive solution for Kk? is 
2 2 2 
k? = em? (1 + Boy + om? Al + P+ Beau? 
2M M 4M P 
vee 2 2 Kk? < 22 
when i 2 ata p-7- > M +m 
p a 
v.28 q7 s 0 


Lori 


abin st x lari 


2c Sno. 


ae 


tue 

Ss 
a 

it | 


a ] is 


Sa 96k -btodiee 


: oe 


=f Uy > 


_ 


sare 


99 


But since case (i) is restricted to k < 2p, and 


since Ms > 2p, one is left with the case Ke < Mot therefore, 
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APPENDIX B: TRANSFORMATION INTO MOMENTUM SPACE 


Rewrite eqns. (IV.14) and (IV.24): 
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The free configuration space wavefunction is [eqns. (IV.16) 


and: €1V.17):1 
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Consider the rhs of (B6); when (B8) is inserted, one is led 


to, apart from the self-energy: 
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ce: «74 25 2 >> > > 
= 27M | d-q oe cd (H(aQ) +0, (K)) 
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which is precisely the term that appears in eqn. (IV.28). 
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APPENDIX C: PARTIAL WAVE DECOMPOSITION 


Define, for the sake of brevity, functions 


V,, (44k) as follows 
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The wave function (or scattering amplitude, apart from a 
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because the p-dependence in V,, (ask) is purely radial. 
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Rewriting eqn. (IV.29) in terms of V. (ark) ies. 
>> >> 
F(p,k) = 20M | d°q Iv, (ark) - vy (rk) I (C6) 
it is obvious that F(p,k) is a function of only the 
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where the following relation has been used: 
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Next, one wants to obtain an equation for the s- 


Waveupart Of then l—matrix, (whien, 1s aefined ‘by (iV .31); 
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This equation therefore serves as the definition of the 
component ~ in terms of ye. Using (C14), one can derive 
an equation for the partial wave T-matrix, qt; and this 


is: 
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Eqn., (C15) ‘holds.for all, even. values of 2; 
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This is the equation that is discussed in Section IV. 
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APPENDIX D: ALGEBRAIC AND NUMERICAL METHODS 


In eqn. (IV.59), the following integral appears: 
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technique as such carries through for the pseudoscalar 
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by Be rn where Me ane «thesvalues of k at the grid 
points j. The left hand side of (D2) is then written as 


apematcrix proauct 


} M, (K5 043) T, (a; oP) 
whereas the right hand side is a column vector SMOSH 12 
The technique then amounts to a matrix inversion, i.e. 
the matrix Si uss Ab is inverted. 

Three terms go into the elements of the matrix Mes 
viz. the three terms on the lhs of (D2). The first two 
are easily obtained, but in order to obtain contributions 
from the third, the integral itself has to be evaluated. 
This is done as follows: 

First, one rewrites it, suppressing all but the 
q-dependence, i.e. 
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Next, the range of integration is divided into N equal 
intervals; i.e. the range is broken into the points 


0 = to tis pecccy wee ty where 
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This interpolation therefore removes T(q) from the 
integral, and hence enables the q-integration to be per- 
formed. Obviously, the method is only good if one assumes 
that T(q) is smooth and singularity-free. Without such an 
assumption, the interpolation is a very bad approximation. 


A typical term in the sum in (D4) is 
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R. (a,b) = | dq war P+ 4-N (D6) 
j Be qtp 


The integrals ee can be evaluated analytically 
as follows: w(q) contains the potential, which in the 


scalar theory is given by (IV.58) 
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In principle, the quartics that appear within the log- 
arithm can be factored, and the logarithm can be replaced 
by sum of logarithms of the types gn| qta| ae 

an| (qtx) 7+y"|. This depends on whether the zeroes of the 
quartic are real or complex. We shall just restrict our- 
selves to the complex case, as the real case is just a 
limittwof ntsc: 


Therefore,.i\(D6)> consists’ of terms ,of isthe ‘type 


b 
a d 2 2 N : 


(D7) can be further simplified by decomposing sty: sO 
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Part of (D8) can be written in the form of Euler's di- 
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by appropriate transformations. s can be real or complex 
depending on whether the zeroes of the quartic are real 
or not. The remainder of the integral (D8) can be done 


in closed form. By this breakdown, therefore, integrals 
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nse may be reduced to combinations of dilogarithms (of 
real and complex arguments) and integrals that can be 
done in closed form. This in turn implies that I. may 
be reduced to a eae form. 

This whole procedure hinges upon being able to 
factor the quartics i.e. determining x,y (or a). This 
is done numerically; the position of the zeroes will of 
course depend upon the value of p. Our region is below 
the single pion threshold. 

Bon.) (DZ). can therefore be written out in matrix 
form. However, the discussion so far has been for k # 0. 
For k = 0, the potential is modified drastically, and the 
integrals no longer are the same. In facts they become 
quite a bit simpler, and the analogues to R. are exactly 
evaluated. No dilogarithms are required in this case. 

The rest of the calculation has to be done numerically. 


Computer programs were written to: 


a) determine the zeroes of the quartic 

b) evaluate the dilogarithm functions (for both real and 
complex arguments) 

c) combine the separate terms 

d) incorporate the k = 0 limit 


e) invert the matrix to determine T. 
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Numerical monet etency was checked by performing 
the matrix inversion with matrices of different sizes. 
Grid sizes of 15 and 20 were picked i.e. the matrices 
were 30 x 30 and 40 x 40. The numbers that are obtained 
are the real and imaginary parts of the half-off-shell 
T-matrix elements and the real and imaginary parts of the 
phase shift, at each value of momentum p (or, in, fact, 
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APPENDIX E: BOUND STATE INTEGRALS 


In Section (IV.6), three integrals were constructed, 
Walz The I, and I.- The first of these can be evaluated 
in closed form. Ie and I, are far more complicated, and 
this appendix outlines the method of handling how The 


integrals are 


00 2 
Se, k“dk 
I, 7A | y zg fy {k) 
OF (kee ) 
Pee 2 peice gi pa cdkwen| emada : 
ea 2 aes 3 Ww 
22m) OS(s oa) 0(q +0) 


(q?+k?+2ME) 7-4m? (q-k) 2-4M?m? 


x kn 
(q7+k?+2ME) 7-4M" (q+k) 7-4m7m* 


With the definition of F, (kK), I_ may be written as 


i) 


F 2 
Pies tatiieee (et et) (El) 
Ss M> 2(2n)* ay Z | 
with 
_ {= kdk - (q?+k7+67) 7-4 (q-k) 7-407 
eae aA riod a yereraeny. we) ee eC hack TREE. 3 Be) 
0 (k° +a") 0 (q° +k°+B8" ) “-4 (qt+k) “-4) 
00 00 26 eZine 2 2 2 
kdk (q°-k“) “-4 (q-k) “-4) 
T - - qdg fn (E3) 
2 Jo (k24+02)7% J0 (Grate Ack) ane 


In arriving at (E2) and (E3), transformation into dimensionless 
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units has been made i.e. 


(q,k) + (Mq,Mk) 


Also, u, B and dX are defined by: 
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For the time being, we shall confine ourselves to I- 
The k-integral in I, (Si say) may be integrated by parts to 
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By interchanging u and v in the second term within 
the integral, it can be combined with the first one, and 


eventually I, can be written in the following form: 
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Consider the u-integral in (E5); because integration goes 
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The structure of the integrand, and also the values 
of the zeroes show that the result is going to be exceedingly 
messy. However, it is obtainable in closed form; then the 
v-integral is evaluated numerically. 
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is slightly different in structure, can be evaluated in an 


I, can be treated identically to I,- Also, Hs which 


analogous fashion. Both I. and Ne were obtained this way, 
and have been tabulated in Tables I to III, for certain 


values of g*. 
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We can go a little further, and examine more closely, 
thereftect of seltoenergysonsthe two-particle system. If 
one looks at the equation for any channel, e.g. the (0,1) 
channel (eqn. (V.43)), the self-energy is seen to be present 
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The zeroes of the first quartic in eqn. (F8) are: 


a ge: w (p) = 2 (M2 43 p*) ar 2 At toMem te 2Mi > 


and therefore, closing the contour in the top half plane, 


only zero to contribute will be 


gq = iw_(p) 


For the second guartic, the zeroes are 


Cree 2a aye Saag +Mm 


Then, (F8) can be integrated, and gives 


27 


* “wt 
a J 


sty 3 


~ ” 

bans 
4 a 
a oF ve | 

{ vy ( ‘ Cre) 
o 4 

1 ~ ube ae 
rp 
£ i 

¢ 
cag 


[i . es 


i@ 


ster oft ots: tsp ‘bnobse 


‘2 
> 


: an 


mOzt ReaneD no tinder nes ytao edt 


bas ,betsxpedni od uD ony ad 


< =] 
ei 
os 


203 
352 a hee 
12M fo (M M aN 


G(p,o) = - eon. ————————s=- Sr 


2 2 
AMe (lie BS) 
M 


Whreh,) tOnorder aae becomes 


2 
or 3 
Ont 5 12 
Cieeews: alam * 2 P 
i.e. 
G(p,o) Es LS <2 m2 m 
p- 4m “oO M 
and, eventually 
1 
4n M fe) 


[3 


This gives the zero-energy renormalisation. 
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Case (b): 


In a similar manner, the infinite limit renormali- 
Sation may be obtained. In this case, however, it can be 
shown that G(p,k) is zero, because the zeroes involved 
are real, and hence do not contribute to the principal 


value integral, i.e. 
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k00 
There is, therefore, no renormalisation in the infinite 
limit. For-large momenta, it is apparent that the approxi- 
mation G(p,k) = 0 is . reasonably good. This approximation 
has been made in obtaining the Born term (Sect. V.4). 

In passing, it should be noted that an analogous 
calculation may be performed for the scalar meson theory, 
and again it can be shown that the zero energy renormali- 
sation is small (for the prescribed value of the coupling 
constant); also it goes to zero in the infinite momentum 
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APPENDIX G: TWO-NUCLEON SCATTERING STATE EQUATIONS 


The field eqn.. (V.3) may be relabelled as follows: 
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APPENDIX H: ISOSPIN DEPENDENCE 


Before one can evaluate the combinations of Clebsch- 
Gordan coefficients that appear in (V.16), one needs the 


following coefficient: 
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Eqns. (H2) to (H5) permit one to explicitly handle the iso- 


spin dependence of (V.16); one has to consider the cases 
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It is obvious from (H8) that the T = 0 state couples 
only to T = 1 states. Also, (H6) and (H8) aa, that the 
lhs and rhs of eqn. (V.16) for T = 0 have the same Clebsch- 
Gordan coefficient. So the equation for the T = 0 channel 


(total spin S being unspecified) is: 


2 2 
Ve. V 
Flys) x Zz Sd ame iah 
<s3v|Ss,> (iss + a ae Toker (Xo, Zy x) 
| > SO > > 
= 3(2T) <s3v|Ss.> farx* We (XK!) Tey ick Zt) 


t v i] v ' Ws t w ier > 
+ 3£5 Og Ss | Axgd x'd"k d“k,d~z' <s}6/S's3><S's5|sy'y 


VATE k 3 3 3 3 
ak 2 3 


239 


ad 


7) Oey 
; ve Tere 
* am aes | te bal 
yf ; “g , ae 7 a Y vent ‘ 
a reae uy i , 
a rt 


ee Pe aa 
11, B<00 ES bo 1 8h Sd MB i ees 


t 


= 
ne 
} coe Oh" : 


the or : : A 7 j 3 
Oa, L eagle, kt 
4 * ; i, — 4 - . a a y 


i ~ A : 
4 
% 
Laat | a 
; ( 
7 t 
i 
a} 7 
ray : y ‘ 4 
tl L, Sol aden al * : 
i é Pisses 5 
ot ' aa" 4 wk ‘ 


a 


7 


8H) Bes (0H) ,@atA ,.aeadeze. I 


- 4 
5 Siste 0 = T oft tedd (8H) moml seozv ft 


’ | y ¢ oO - ~ | a pu ; - 2 
) omspe sddt Sven O = T rot (eL.V) ape 20 aaa @ 
lias 


yon 


ote a) io ae Sr ste: tt. bee eens eel ~ #2 
‘fT! y T sdt ~wol monzagps Ssht on 


\ 
;2ik (belitioesqeny pre 
u¢ 
~ 
f 4 
'' *¥ mea . 
- ~ oe — 9 7 
{ Bes \ ns (—-- ae + —wwt ) < ae 
‘e. oe mS £ 
cr i] = ame 


*« =< es = tet ae ey es 
O*xt ak 4 ~ ('x—) Wok 8S im 7 
: © ed ais | | Se - 
so 


aeee tees be ne ~ 18 sy bie —— ‘ Mts A a a . . . ' . 
i" ~at._2 c,- * pe ¢ Cs i ve j a Oo mw fe} at 5 x w & i 
: a - . od i ¢ : b : A on 
e% ~ 
= : 


140 


S 
, we 
x <s} 6|Ss,>F, 


+ the corresponding bound state term. (H9) 
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An analogous treatment to that in (a) can be done 
here. Since the method is identical, and the algebra much 
longer than in (a), it will not be reproduced here. Instead, 


the final equation will just be written down, and is: 
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+ the corresponding bound state terms (H10) 


Unlike the T = 0 case, the T = 1 state couples to 
both T = 0 and T = 1 intermediate states, with different 


weights. 
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APPENDEX 41: SPUN »DERENDENGE 


The calculation here is straightforward, and is 
very similar to that in App. H. A relation identical to 
eqn. (H2) is obtained for the spins. The full expression 


required is 
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Therefore, as before, (11) becomes, for S = 0: 


' * 
Sw) (eae) ae s 8 14t8| 00>F ge 


Or 


’ '* ' Ix 
i cits Glee mana Dy. eee (leprae saat 
aie om F 


fe 
LS) 


OT 


1 tx 
sa) CaDNR pk ohn er [5 5 51476|002F 0 pe 
§ 


5 


142 


(iT) 
ee 
ae CS 


a’ 


WITH 


: T) 


“Hi i 


7 
tle a 


yee ‘* Vwi bes 
T "By 


wo ad grey 


siom sit we 


q 


2SmooSd 


.. « 
00} d~>4 45 
oat 
$b 


ae 
sOREGR, pi: 


t 


ii | 


vd bes iqox ois. 
he f, 


(£T) 


sec weg 


re 


\Stoted 46 \s 


UX 
" - S244 <-6|00>F27 Ul OT 
3 
oa (BF), GT) _4(-5 oan <+6| 00> Bent Ode 
5 3° 
q 1 oe 
+ 5 Sot <+6|00>Ft7 rit mee 
3 
| t *& 
i 6-6-5 5 §1_<-5|00>F0" FOF 
5 35 
' tx ' '* 
peer ye8 20 G0>Pis Pe 64 Berl 00M SES 
2 sie Si+ OT 
3 3 F 
(I2) 


In order to evaluate (I2) one needs the following: 
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This is the term corresponding to (16); once again, as in 
the case of isospin, it is apparent from (16) and (I7) that 
singlet states (S = 0) couple only to triplet states (S = 1) 


while triplet states themselves couple to both. 
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APPENDIX J: ETACR AND COMPLETENESS 


Eqns. (V.5) for the ETACR leads to a completeness 


relation for two-nucleon states in the following manner: 
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